Anisotropic properties of spin avalanches in crystals of nanomagnets 
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Anisotropy effects for spin avalanches in crystals of nanomagnets are studied theoretically with the 
external magnetic field applied at an arbitrary angle to the easy axis. Starting with the Hamiltonian 
for a single nanomagnet in the crystal, the two essential quantities characterizing spin avalanches are 
calculated: the activation energy and the Zeeman energy. The calculation is performed numerically 
for the wide range of angles and analytical formulas are derived within the limit of small angles. 
The anisotropic properties of a single nanomagnet lead to anisotropic behavior of the magnetic 
deflagration speed. Modifications of the magnetic deflagration speed are investigated for different 
angles between the external magnetic field and the easy axis of the crystals. Anisotropic properties 
of magnetic detonation are also studied, which concern, first of all, temperature behind the leading 
shock and the characteristic time of spin switching in the detonation. 
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I. INTRODUCTION 

Single-molecule magnets (nanomagnets) embedded in 
crystals are compounds that exhibit unique physical 
properties with promising applications to quantum com- 
puting and data storageii"— In particular, they can pos- 
sess a large effective spin number {S ^ 10)^ and show an 
anisotropy with respect to the orientation of this spin, 
with the lowest energy corresponding to an "easy axis" 
of the crystal)^ That is, the potential energy as a func- 
tion of the orientation of the spin exhibits a double- well 
structure, even in the absence of any external magnetic 
field. In the presence of an external magnetic field par- 
allel to the easy axis, the two wells are asymmetric and 
the spin aligns with the field. Upon a sudden reversal 
of the field, the internal crystal anisotropy creates a bar- 
rier to the flip of the spin, and relaxation may take place 
through spin tunneling^— It is also possible to trigger 
locally the relaxation and, as it releases energy, observe 
the propagation of a spin reversal front, corresponding to 
a magnetic deflagration!^— or detonationr^'^ depend- 
ing on the speed and structure of the front. Magnetic de- 
flagration and detonation have much in common with the 
respective combustion phenomen a^'^i^'^ (including the ter- 
minology) ; there are even indications on the possibility of 
magnetic deflagration-to-detonation transition similar to 
that studied intensively within combustion science.—"— 

Up to now, the research on magnetic deflagration and 
detonation has mostly been restricted to unidimensional 
models, where the external magnetic field is co-linear 
with the easy axis and the spin avalanche front prop- 
agates along the same axis. Within such a restriction, 
one obviously loses the possibility of an anisotropic spin 
interaction with the magnetic field, together with an 
anisotropic propagation of the avalanche fronts. Al- 
though the importance of and interest in the anisotropic 
properties of spin avalanches was expressed from the 
very beginning;^ only a few papers addressed these 
properties which may be explained by the experi- 



mental difficulties encountered in its study. In partic- 
ular, Ref. [2^ investigated experimentally the possibility 
of spin- avalanche initiation ( "ignition" ) for the magnetic 
field inclined at an arbitrary angle to the easy axis. In 
Ref. ^Sd, the authors compared the magnetic defiagration 
speed for propagation along the easy axis (c) and the hard 
axes (a or b) with the magnetic field collinear with the 
front velocity vector. Thus, although the experimental 
data on the subject is limited, the anisotropic properties 
of the magnetic deflagration and detonation may be in- 
vestigated using nanomagnet model Hamiltonians»i^ To 
the best of our knowledge, no theoretical investigation of 
these anisotropic properties has been performed so far. 

In the present paper, we explore the effects of mis- 
alignment between the external magnetic field and the 
easy axis. We shall focus on the development of a model 
for magnetic deflagration and detonation in a crystal 
of single-molecule magnets in a generic magnetic field. 
While this model can be applied to any such system, spe- 
cific calculations will be based on Mni2-acetate, which 
has an effective spin number S = lOJ^i^ Starting with 
the Hamiltonian for a single magnet embedded in the 
crystal, we calculate the two essential quantities - the ac- 
tivation energy and the Zeeman energy - characterizing 
the spin avalanche. We investigate modifications of the 
magnetic defiagration speed produced by misalignment 
of the magnetic field with the easy axis. We also study 
the anisotropic properties of magnetic detonation, focus- 
ing on the temperature behind the leading shock and 
for completed spin reversal, and the characteristic time 
of spin switching. Unlike for magnetic deflagration, the 
magnetic detonation speed is determined by the sound 
speed and does not depend on the direction of the exter- 
nal magnetic fleld. 

The paper is organized as follows. We start by present- 
ing, in the next section, the quantum-mechanical calcula- 
tion of the activation energy and the Zeeman energy. We 
then derive, in Sec. IIIIl approximate analytical formu- 
las for these values, based either on quantum-mechanical 
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perturbation theory or on a classical model for the spin. 
In Sec. IIV[ we consider the implications of the quantum- 
mechanical results on magnetic deflagration and deto- 
nation properties. Finally, we summarize our results in 
Sec. El 



II. QUANTUM-MECHANICAL DERIVATION 
OF THE ACTIVATION AND ZEEMAN 
ENERGIES 

A. Hamiltonian for a single-molecule magnet 

A rather elaborate spin Hamiltonian for a molecular 
magnet, such as Mni2-acetate, can be written as^ 



n = -DS'l - BSt 



9f^B 



Ht (^cos (f>Sx + sin 

4 , ^4 



s^) + c (5+ 



n', 



(1) 



with the spin raising and lowering operators S± = 
iSy. The first two terms of Eq. ^ correspond to the 
uniaxial magnetic anisotropy, while the third term is the 
interaction with a magnetic field H, oriented along the 
spherical angles (0,0), with the components 



H sin 6 cos ( 



H cos 9, 



while 



(2) 



(3) 



is the transverse magnetic field. The 4th and 5th terms 
of Eq. ([ij are transverse anisotropy terms (inherent to 
the molecule), and H' contains additional terms due to 
the inter-molecular dipole interaction and the hyperfine 
interaction with the spin of the nuclei. A set of values 
for the parameters in this Hamiltonian for Mni2-acetate 
can be found in Tab. HI 

Even in the absence of a magnetic field, the presence 
of the transverse anisotropy terms makes it such that the 
eigenstates of Sz are not eigenstates of the full Hamilto- 
nian ((J). Nevertheless, due to the small values of E and 



TABLE L Values of the different parameters of the spin 
Hamiltonian, Eq. for Mni2-acetate. 



Parameter 



Value 



Ref. 



g 


1.93 




D 


0.548 K 


[16] 


B 


1.17 X 10"^ K 




E 


1.0 X 10"^ K 


fl6] 


C 


2.2 X lO^"^ K 


[16] 



C, it is still informative to discuss the problem in terms 
of the magnetic quantum number associated with Sz ■ 
We plot in Fig. [TJa) the energy of the Mz eigenstates in 
a field of F = 1 T aligned along z (6* = </> = 0). An E 
and C are small perturbations, the eigenvalues of % are 
almost those of Sz, and only the Mz = 10 level is signif- 
icantly present in the ground state. Rotating the polar 
angle to 6* = 7r/3 changes not only the energy of the Mz 
levels. Fig. [TJb), but also increases the "population" of 
the different Mz in the ground state of the system, that 
is, the projection of the ground state ipg on the eigen- 
states of Sz, \{Mz I V'g)!^- While the ground state is still 
located close to the maximum projection of the spin on 
the z axis, i.e., the system is found in a single well of the 
double-well structure, the energy of the ground state is 
higher than that of the lowest Mz level. This can be ob- 
served by considering the expectation value of {Sz) in the 
ground state of Hamiltonian ^ for different orientations 
of the magnetic field. Fig. [2l 

The combination of the change of the level structure 
and the projection of the initial and ground states on 
many levels will affect the values of the activation and 
Zeeman energies, as described in Sec. IH Bl In all cases, 
we need the anisotropy to play the dominant role, so that 
the double-well structure of the spin energy is present. 
Defining the anisotropy field as^i 



Ha = {2S - 1) 



D + B[2S[S-1) + 1] 



91J-B 



(4) 



we must have at all times Hz < Ha and Ht ^ Ha, with 
Ha ~ 11.1 T for Mni2-acetate. Also, while the Hamil- 
tonian ^ is different along x and y, this leads only to 
minimal modifications in the energy as (j) is varied, and 
we will thus concentrate on the behavior in the xz-plane, 
i.e., for (h — 0. 



B. Determining the activation and Zeeman 
energies 

The physical situation we consider here is the follow- 
ing. Initially, a crystal of molecular magnets is immersed 
in an external magnetic field H_ , which is then very 
rapidly inverted to a new field H — — H_. Because of the 
magnetic anisotropy, the system is then in a metastable 
state, and an energy barrier must be overcome for it to 
relax to the new ground state. The relaxation of a given 
molecular magnet can then happen through spin tun- 
neling, where less energy than the barrier height is re- 
quired, or by thermal excitation above the barrier. The 
molecular magnet thereby releases the thermal energy 
equivalent to the difference in energy between the ini- 
tial metastable state and the actual ground state. This 
thermal energy can then contribute to the relaxation of 
neighbouring molecular magnets, hence the possibility of 
deflagration and detonation inside the crystal. 

In order to serve for the study of deflagration and det- 
onation, our model must therefore produce two main val- 
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FIG. 1. (Color online) Energies of the eigenstates of Sz, labelled by the quantum number Mz, in an external field of 1 T 
oriented along (a) = (fi = 0; (b) 6 = tt/S, (f) = 0. The projection of the ground state on the different Mz levels, |(Afz | ^g)!^, is 
schematically represented by the thickness of the line (in a logarithmic-like scale), with dotted lines 0. 




FIG. 2. (Color online) Expectation value of {Sz) in the ground 
state \tpg) of Hamiltonian ([1} for Mni2-acetate as a function of 
the polar angle 9 between the magnetic field H and the easy 
anisotropy ajcis of the crystal, with <?!> = 0, for three different 
magnitudes of the field. 

ues, the activation energy £a, i-S-, the difference between 
the maximum energy of the molecular magnet in the field 
H and the energy of the initial metastable state, and the 
Zeeman energy Q, corresponding to the difference be- 
tween the metastable state and the ground state in the 
field H. Therefore, we first solve 

?i„ IV.) = £-,0 \A) , (5) 

with "H- the Hamiltonian using the field H_ (i.e., the 
field before inversion), for £__o the lowest eigenvalue of 
fi-, and then calculate the energy of that state in the 
field H, 

£^ = mn\^^). (6) 



To get the barrier height, we consider the spin-phonon 
coupling as a sum over products of all the spin operators 
Sx, Sy, and Sz, (see, e.g., Ref. 'sT'), such that the system 
overcomes the barrier by stepping through intermediate 
states up to the state of highest energy £niax in the field 
i.e., 

T-C l-^max) = fmax IV'max) , (7) 

such that 

£a ^max £i- (8) 

Note that this model takes into account the effect of tun- 
neling on the position of the energy levels, but not the 
dynamical effects of tunneling. In other words, we con- 
sider that the crossing of the barrier due to thermal ex- 
citation will be much faster than the tunneling across it 
(opposite to what is studied in Ref. [29l). 

The Zeeman energy is itself found from the state of 
lowest energy £mi„ in the field H, 

H IV'min) = £min IV'min) , (9) 

as 

Q = - £,„i„. (10) 

Both £a and Q are easily calculated numerically, and 
some results for a magnetic field in the xz-plane are pre- 
sented in Fig. [3l From the structure of the Hamilto- 
nian ID), while it is clear that these values are mirrored 
about 9 = and 9 = 7r/2, there is a difference in be- 
haviour of the curves around these two angles. While the 
Hamiltonian is symmetric about both 9 = and 7r/2, the 
presence of Ht, Eq. makes the first derivative of the 
energy discontinuous at 9 = tt/2, and this is reflected in 
both £a and Q, as can be seen in Fig. 
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FIG. 3. (Color online) Quantum-mechanical calculation of (a) the activation energy and (b) the Zeeman energy of Mni2-acetate 
as a function of the polar angle 8 between the magnetic field H and the easy anisotropy axis of the crystal, with (f} — 0, for 
three different magnitudes of the field. 
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FIG. 4. (Color online) Quantum-mechanical calculation of the activation energy (solid line) and the Zeeman energy (dashed 
line) of Mni2-acetate around the symmetry angles (a) 6 = and (b) 6 = n/2, for (j) = and H — IT. 



C. Range of validity of the model 

An underlying assumption of this model is that, ini- 
tially, a single quantum level of the molecular magnet 
is populated. This is of course dependent on the initial 
temperature of the system, so it is useful to also look at 
the difference in energy between the lowest state in field 
H_ and the next-to- lowest. We denote this quantity by 
Sgap and a plot of its value can be found in Fig. [5j The 
curves clearly show a change of behavior for a certain 
value of the angle 9, which can be easily understood as 
follows. If, for the sake of the explanation, we neglect 
the fact that more than one level is populated and 
only think in terms of the energies of the states, for 
small angles the difference in energy corresponds to that 
between = — 10 and = —9 (for H_, the structure 
is reversed with respect to figure Fig.[T|). Above a certain 



value of 9, the component of the magnetic field along z, 
Hz, is too weak, such that the level = 10 is actually 
lower in energy than = —9, and Sg^p corresponds 
to the difference between the ortho- and paramagnetic 
states of the crystal. The kink in Sg^p is therefore due to 
the shift from a structure of the type of Fig. [IJa) to that 
of Fig. mb). 



In the first case, where the energy gap is between two 
eigenstates on the same side of the well, thermal excita- 
tion will lead to a small correction of the activation and 
Zeeman energies, as the initial state of the system will 
have a higher energy than calculated here. In the latter 
case, the thermal energy will lead to an initial projection 
on the levels in both wells, leading to a breakdown of the 
model. 
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FIG. 5. (Color online) Energy gap between the lowest and 
next-to-lowest eigenstate of the Hamiltonian ^ for Mni2- 
acetate as a function of the polar angle 9 between the mag- 
netic field H_ and the easy anisotropy axis of the crystal, 
with (j) — 0, for three different magnitudes of the field. 



III. APPROXIMATE FORMULAS FOR THE 
ACTIVATION AND ZEEMAN ENERGIES 

While an implementation of the rescription of Sec. Ill Bl 
rehes on the numerical solution of an eigenvalue system, 
this can be done in real time when coupled to a simulation 
of deflagration or detonation. However, it is also useful 
to have analytical formulas, which can give insight into 
the physics governing the processes. We therefore derive 
approximate equations for £a and Q, for the case where 
the external magnetic field is nearly aligned with the easy 
axis of the crystal, i.e., Ht ^ H^. For this purpose, we 



will also consider the simplified Hamiltonian 

n^-DS^.-BSt-g^iB^Hj^ + HrSj^ (11) 

where we have set cf) ~ and neglected the transverse 
anisotropy terms. 

A. Perturbative approach 

With the exception of the term in Sx , the Hamiltonian 
in Eq. ([TT]) is considered by many authors as the "unper- 
turbed" Hamiltonian, the other terms being responsible 
for a slight shift of the energy levels and for magnetic 
tunneling. This is also the case for Ht <C i?^, so let us 
define the unperturbed Hamiltonian as 

Ho = -DSl - BSt ~ g^isHj, , (12) 

with the perturbation 

n' = -cjiiBHrk ■ (13) 

The eigenstates of the unperturbed Hamiltonian will be 
written as 

Tio |M,) = |M,) , (14) 

with the unperturbed energy 

4? = -DMl - BMt - gfiBH.A'h. (15) 

There is no first order correction to the energy, i.e., 

4}^ = (A4|?i'|Af,) =0, (16) 

since the diagonal elements of Sx are 0. We thus need to 
consider second-order corrections. 



c(2) 
■'M 



E 



\{M'jn'\M,) 



c(0) _ c(0) 



{M, ~ l\n'\M,) 



(0) 
M 



c(o) _ AO) 



V 2 



5(S+l)-M,(A/, + l) 



(2Af, + 1) [-D - B (2Af| + 2Af, + 1)] - g^iBH, 
S{S + l)-M^{M^-l) 



(2A4 -l)[D + B (2Af2 _ 2M^ + 1)] + gfisH^ 
such that the total energy is given to second order by 



_ c(0) 



£m = S 



M 



E 



(1) 

M 



s 



-DMl - BMt ~ gtiBH^Ah 



V 2 



S{S+1)^M,{M, + 1) 



(17) 



(2A4 + 1) [-D - B (2Af| + 2A/, + 1)] - g^ieH, 



+ 



S{S +1) - M,{M^ -1) 



{2M, -l)[D + B (2Af| - 2A^, + 1)] + qiibH, 



(18) 
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Explicitly, we have the energy of the initial state \—S), 



e.s = -nS' - BS' + ,,sB.S + ^ _ ,^ ^ l^^^"";- ^ _ ^^^^^^ . (19) 

and of the ground state | S) , 

after inversion of the field. The Zeeman energy is thus found to be 

Q = £^s-Ss=9t^BH,sh + 2 [m^BHr) ^1 ^21) 

1 [2S-lf[D + B{S^-2S + l)f~{giiBH.f\ 



,2 



Calculating the activation energy is more tricky, as it 
requires knowledge of the value of Mz for which the en- 
ergy is maximum. We remedy this by considering Mz 
to be real, and not limited to integer values. Using the 
unperturbed energy, Eq. p3|) . we find 



maxf 



(0) 
M 



D 



(7/9) 



1/3 



(37) 



1/3 



2B 



where we have defined 



7 = 9B^gfiBHz - 
We also get that 



8531)3 + 2754 {guBHz 



lim A/ir 



2D 



(22) 



(23) 



(24) 



We finally can get the approximate activation energy by 
substituting A/max into Eq. (jl8p and subtracting the en- 
ergy of the initial state [Eq. (|19p ]. i.e., 



(25) 



In Fig. ini we present the relative error on the calcu- 
lation of £a and Q, as compared to the exact quantum- 
mechanical calculation, as presented in Sec. Ill Bl but for 
the Hamiltonian (|lip . As expected, the results are in 
good agreement for small angles, but a strong deviation 
is observed as Ht becomes non-negligible compared to 
i?z- A much better approximation is obtained for the 
Zeeman energy, in great part because levels close to the 
top of the barrier are more affected than those at the 
bottom of the wells and because of the additional ap- 
proximation that is need to determine the value of 
for which the energy is maximum. 



B. Classical approach 

Following the approach of Macia et al.y^ we shall now 
treat the spin of the nanomagnet as a classical vector S. 




FIG. 6. (Color online) Relative error on the activation (A£^a) 
and Zeeman (AQ) energies for Mni2-acetate calculated using 
perturbation theory for the Hamiltonian Eq. for two 

different magnitudes of the field. 



By deriving the dependence of the energies with respect 
to the orientation of the spin vector, it will be possible 
determine the activation and Zeeman energies, following 
the same method as prescribed above for the quantum 
Hamiltonian. 

From the Hamiltonian ([TT|) , we get the classical formu- 
lation of the energy 



.class ^ _^^2 _ ^^4 „ (^^^^ ^ 

= —D {S cos a)^ — B {S cos a)^ 

— gUB {HzS cos a + HtS sin a) , (26) 



where a is the angle between the spin vector and the 
z axis. The minimum energy, from which we can get 
the orientation of the initial spin vector and that of its 
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ground state, is therefore found by solving 



da 



2DS^ cos a sin a + 4_BS"* cos'^ a sin a 
— gUsS {—Hasina + Ht cosa) — 0. (27) 



Making the assumption that the transverse field Ht is 
small compared to the internal anisotropy, see Eq. (|4]), 
we get that the spin vector will be nearly aligned with 
the easy axis (z), and the external field will introduce 
only a slight deviation. This is indeed what is observed 
for the full-quantum calculation in Fig. [2j The angle a 
is thus small, such that we can approximate Eq. (j27p by 



d£' 



class 



da 

and we get 



2DS^a + 4BS^a - g^lBS {-H^a + Ht) = 0, 

(28) 



2DS + ABS^ + g^lBH, 
Thus, the energy of the ground state, £[ 



class 
9 ' 



(29) 



is obtained 
I 



from Eq. ((26)) using Eq. ([29)1 for a. 

The energy of the initial state, £^^^^^ is also obtained 
from Eq. (|26|) . but using the angle of the spin vector 
amin,- in the inverted field, H_. Following the above 
procedure, we easily find that 



an 



TT + q;„ 



(30) 



[The symmetry of Eq. (j26p with respect to the inversion 
of the external field can also be used to demonstrate the 
relation between amin,- and amin-] We finally get 



£. 



:lass 



- £'^ 

= i (a = vr + amin) - £ (" = amin) 

= 2giiBS (i?^ cos amin + Ht sin a^in) • (31) 



To calculate the activation energy, we again need to 
determine the highest energy the spin vector will have to 
overcome as its angle goes from amin,- to amin- Plotting 
Eq. (1^51) as a function of a, one can easily see that the 
maximum is around a ~ 37r/2. Making the substitution 
ckmax = 37r/2 + e, with e a small angle, into Eq. (P7)) . we 
have 



2L>S'^ sin e cos e - 4BS"^sin3e cos e-5^BS'(i?2 cos e + iJr sine) w ~2DS'^e- ^^BS'^e^ ~ g^isS {H^ + Htc) = 0. (32) 



Solving for e, we get 



37r 1 



2 251/35 



-gfisH^ 



V 



i2DS + gfiBHT) 



{2DS + gfiBHrf 
27B53 



27BS^ 



(gt^BH^) 




From this expression for amax, we find the corresponding energy £^1^.^^ leading to 

cclass cclass (rclass (rclass 



cclass cclass 



%a 

„2 



-^max 
rclass / 



nax) - £'' '''"'(a = TT + amin) 
DS^ (COS^ amax " COS^ amin) " B (cOS"* amax ~ COS"* amin) 
gfiBS [H^ (cos amax + COS amin) + Ht (sin amax + sin amin)] 



(33) 



(34) 



We once more calculate the relative error with respect 
to the exact quantum-mechanical values using the Hamil- 
tonian ([TT]) . see Fig. [T] The result is markedly better 
than that obtained using perturbation theory (Fig. [SJ, 
even for greater values of the angle 9. This can be easily 
explained by the fact that the classical model is based 
on a much different approximation, namely that the spin 
only slightly deviates from being aligned with the easy 
(z) axis. This gives a validity over a much greater range 
than what was given by treating the transverse field Ht 
as a perturbation. 



IV. ANISOTROPIC PROPERTIES OF 
MAGNETIC DEFLAGRATION AND 
DETONATION 

A. Magnetic deflagration 

In this subsection, we investigate the magnetic defla- 
gration speed for an arbitrary angle between the mag- 
netic fleld and the easy axis; the next subsection will 
be devoted to magnetic detonation. We stress here that 
the propagation of magnetic deflagration involves four 
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FIG. 7. (Color online) Relative error on the activation {A£a) 
and Zeeman (AQ) energies for Mni2-acetate calculated in the 
classical approximation, based on the Hamiltonian Eq. (|lip . 
for two different magnitudes of the field. 





important vector values: the magnetic field intensity H, 
the magnetization M, the temperature gradient VT, and 
the heat flux kVT (with k being the tensor of thermal 
conduction). The latter two are in general not parallel 
because of the crystal anisotropy. Among these values, 
the temperature gradient VT determines the direction of 
front propagation, while the magnetic field intensity H 
and magnetization M specify the activation energy and 
the Zeeman energy of the spin reversal as discussed in the 
previous sections. We stress that the vectors H and M 
are not related to the direction of the deflagration front 
velocity, but infiuence the absolute value of that velocity. 
We also point out that 1) the direction of the magnetic 
field H is controlled by the experimental set-up; 2) the 
direction of the magnetization M correlates strongly with 
the easy axis (c-axis) of the crystal (see the calculations 
above and Fig. [5]); 3) the direction of the temperature 
gradient VT and front propagation are determined by 
the ignition conditions, e.g., by surface acoustic waves;?--^ 
and 4) the direction of the heat fiux k\7T results from the 
anisotropic thermal conduction of the crystal. The differ- 
ent directions defined by these four vectors open a wide 
parameter space for experimental studies of anisotropic 
crystal properties, both magnetic and thermal. As an 
example. Fig. |8] illustrates some possibilities of the mag- 
netic defiagration geometry [because of the small factor 
E = 10~^ K in the Hamiltonian for Mni2-acetate (see 
Tab. m the difference between the a and b crystal axes 
is minor] . Figure EKa) shows the commonly investigated 
case of a defiagration front propagating along the easy 
axis with the magnetic field and magnetization aligned 
along the same axis. In Fig. |Hl[b), the magnetic field 
points along the easy axis, but the magnetic defiagration 
front propagates along the hard axis (axis a or b). Ob- 
viously, both the activation and Zeeman energies are the 
same for the geometries of Fig. (SJa) and (b); but the de- 



(c) deflagration front 

FIG. 8. (Color online) Schematic of the deflagration front ge- 
ometry in the crystal of nanomagnets for the following three 
cases: (a) the external magnetic field and the front propa- 
gation are parallel to the easy axis; (b) the magnetic field 
is parallel, but the front propagation is perpendicular to the 
easy axis; (c) front propagation is parallel to the easy axis, 
but the field is at some angle to the axis. 



flagration speed is different because of different thermal 
conduction along the easy and hard axes. In particu- 
lar, by comparing the magnetic deflagration speed for 
these two geometries, C^/(a,b) and C//(c), one can measure 
the ratio of the thermal conduction coefficients Ha,b/f^c 
quantitatively as Ka^bji^c — [C^/(a,b)/t^/(c)]^; a qualita- 
tive comparison of these coefficients was performed in 
Ref. H^- Finally, Fig. |5Jc) shows the geometry with the 
front propagating along the easy axis, but with the mag- 
netic field directed at some arbitrary angle to the axis. 
In this section we focus on the geometry of Fig.lSJc). For 
large magnetization values, M ^ H , this geometry in- 



9 



volves refraction of the magnetic field at the deflagration 
front. Still, for the crystals of nanomagnets used in the 
experimental studies so far, the magnetization is small, 
M ^ H, and the refraction effects may be neglected. In 
principle, one may consider an even more general geome- 
try than that shown in Fig. [Sjc) with both the magnetic 
field and the front speed aligned at some angle to the 
easy axis. However, at present there is no quantitative 
experimental data for the ratio Ka,b/i^c', therefore, such a 
general case involves unidentified parameters and, with- 
out proper experimental support, it may be considered 
only as an hypothetical study. 

Within the geometry of Fig. ^c) , the governing equa- 
tions for magnetic deflagration arei^>2S 



'dt 



and 



dn 
'dt 



TR 



■ exp 



V ■ (kVS) ~ Q 



' T 



dn 
'dt 



1 



cxp (Q/T) + 1 



(35) 



(36) 



where £ is the phonon energy, T is temperature, n is 
the fraction of molecules in the metastable state (i.e., 
normalized concentration), tr is the coefficient of time 
dimension characterizing the kinetics of the spin switch- 
ing. We also take into account here the possibility of 
a non-zero final fraction of molecules in the metastable 
state in the case of relatively low heating (low Zeeman 
energy), which has been termed "incomplete magnetic 
burning" in Ref. [l^. This fraction is given b y^^'^° 



1 



exp (Q/T) + 1 ' 



(37) 



which is (obviously) taken into account in Eq. p6|) : here 
the label / refers to the final state of the system after 
the avalanche. As we can see from Figs. [3] and [TUl the 
concentration Uf cannot be neglected in the case of a 
small magnetic field and/or strong misalignment with the 
c-axis. The phonon energy and crystal temperature in 
Eqs. (l35t and (|36)) are related according to ^^i^^ 



Akp 



a 



1 



T 



(38) 



where A = 127r^/5 corresponds to the simple crystal 
model, ks is the Boltzmann constant, a is the problem 
dimension (we take a = 3, as we consider the 3D case), 
Qd is the Debye temperature, with B/j = 38 K for Mni2 
acetate. The thermal conduction may also depend on 
temperature; Refs. [Tol and[20l considered the dependence 
in the form n oc with the parameter (3 within the 
range to 13/3. Below we show that the case of con- 
stant thermal conduction, i.e., /3 — 0, gives the best fit 
to the experimental data.-i^ 

We consider the stationary solution to Eqs. (1551) and 
(I36p for a planar magnetic deflagration front propagating 
with constant velocity Uf along the z-axis (the easy axis). 



In the reference frame of the front, 
reduce to 



Eqs. dsn and ^ 



d£_ 

dz \ dz 



Uf- 



dn 



TR 



■ exp 



£a 

' T 



exp {Q/T) + 1 



(39) 



(40) 



The boundary conditions for the system are determined 
by the initial energy fo (temperature To) far ahead of 
the front, and the final energy £f (temperature Tf) far 
behind the front. The initial and final energies (tempera- 
tures) are related by the condition of energy conservation 
£o + Qna = £f + Quf, or 



at; 



(« + i)0S 



Q 1- 



1 



AT^^ 



exp (g/To) + 

Q 

1' 



-, (41) 

(a + 1)9^ exp(g/r/) + l' ^ ^ 

which follows from Eq. pop . In particular, our calcula- 
tions use a low initial temperature. To = 0.2 K, which 
allows reducing Eq. (PT|) to the simpler form 



at; 



Q/Tf 



(a + 1)0'^ 1 + exp {-Q/Tf) 



(42) 



We calculate final temperature Tf and the final molecule 
fraction in the metastable state n / numerically for differ- 
ent strengths and inclinations of the magnetic field; the 
results obtained are presented in Figs. and [TUl together 
with the scaled activation energy £a/Tf, which plays an 
important role for the deflagration front dynamics. As 
we can see, the temperature Tf increases with the held 
and decreases with the angle; the scaled activation en- 
ergy £a/Tf decreases with the field and increases with 
the angle. Still, this decrease/increase is not dramatic; 
for example, for iJ = 1 T, the temperature changes from 
12.2 K to 6.0 K and the scaled activation energy from 4.5 
to 11.8 as the angle 9 varies from to 7r/2. We will see 
below that the variations of the deflagration speed are 
much stronger because the speed is sensitive to both the 
final temperature and the scaled activation energy. 

A qualitative understanding of the magnetic defla- 
gration speed may be obtained from the Zeldovich- 
Frank-Kamenetsky theory, from which we have the 



expression 



19.33 



i^f 
Ztr 



exp 



£a_ 

2Tf 



where Z is the Zeldovich number. 



Z = 



£a Q{l- nf) 
Tf CfTf 



£a 



(a + l)T/ 



(43) 



(44) 



and Cf = {d£/dT)f is the heat capacity in the heated 
crystal. The flnal relation in Eq. becomes an accu- 
rate equality for the case of complete magnetic burning. 
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FIG. 9. (Color online) The parameters of the magnetic de- 
flagration vs the angle between the magnetic field and the 
easy axis: (a) final temperature, (b) final concentration of 
the metastable molecules, (c) scaled activation energy. 



n/ 1. The Zeldovich-Frank-Kamenetsky theory, giv- 
ing the speed [Eq. (|43)) ]. holds only for large values of 
the Zeldovich number Z ^ 1. Such large values are com- 
mon in combustion problems but rather unusual for 
magnetic deflagration. As we can see from Figs. 1^ and [TUl 
the Zeldovich-Frank-Kamenetsky theory may be applied 
to magnetic deflagration only for the cases of sufficiently 
low field and high angles between the magnetic field and 
the easy axis approaching 7r/2. In the case of a moder- 
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FIG. 10. (Color online) The parameters of the magnetic de- 
flagration vs the magnetic field: (a) final temperature, (b) 
final concentration of the metastable molecules, (c) scaled ac- 
tivation energy. 



ate Zeldovich number, as often encountered in magnetic 
deflagration, the deflagration speed may be calculated 
numerically on the basis of Eqs. (1391) and (|40)) using the 
numerical method of Refs. and |34 

We point out that the problem contains a number of 
parameters, and still remains a challenging task for ex- 
perimental measurements, such as the thermal conduc- 
tion Kf and the coefficient of time dimension character- 
izing spin-switching tr. The temperature dependence 
of thermal conduction n oc is also unclear, with the 
factor /3 treated as a free parameter in Refs. [l^ and [s^ 
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FIG. 11. (Color online) Comparison of the experiments and 
numerical calculations for the magnetic deflagration speed 
versus the applied magnetic field. The markers show the ex- 
perimental data of Ref. [itI . The solid lines present the nu- 
merical solutions for different temperature dependencies of 
the thermal conduction coefficient with /3 = 0; 3; 13/3 and 
K-f/fR = 207m^/s^ providing the best fit of the experimental 
data. The dashed line presents the analytical theory, Eq. 
plotted for /3 = and Kf/rR^ 207 mVs^. 
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changing within the range ofO</3<13/3. We suggest 
here choosing particular values of the unknown parame- 
ters by comparing numerical results to the experimental 
dataii^ obtained for the magnetic field aligned along the 
easy axis. Figure [TT] presents the magnetic deflagration 
speed versus the magnetic field calculated for different 
values of Hf/TR and /3. Comparison to the experimental 
data suggests the parameter values i^f/rR = 207 m^/s^ 
and /3 = 0, which provide the best fit for the experimen- 
tal results (red line) and which we use in the following 
for investigating the anisotropic properties of magnetic 
deflagration. The method of least squares was used to 
fit the data. As we can see in Fig. [Til a strong tem- 
perature dependence of the thermal conduction n cx 
with /3 = 3, 13/3 leads to an excessively strong depen- 
dence of the deflagration speed on the magnetic field, 
which does not reproduce the experimental data prop- 
erly. Figure [TT] shows also the analytical predictions of 
the Zeldovich-Frank-Kamenetsky theory, Eq. (jTS)) . plot- 
ted by the dashed line for the same parameters k / /r^j = 
207m^/s^ and /3 = as the numerical solution. As we 
can see, the analytical theory provides only qualitative 
predictions in the experimentally interesting parameter 
range. 

The numerical results for the magnetic deflagration 
speed are presented in Fig. [T2j (a) versus the angle be- 
tween the magnetic field and the easy axis for different 
strength of the magnetic field; (b) versus the magnetic 
field strength for different values of the angle. All plots in 
Fig. [T2| demonstrate the same tendencies - monotonic in- 
crease of the deflagration speed with the field and strong 
decrease with the angle. The tendencies are qualitatively 
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FIG. 12. (Color online) Magnetic deflagration speed (a) 
versus the angle between the magnetic field and the easy axis, 
and (b) versus the magnetic field strength for different angles. 



the same as one had for the final temperature; still, they 
are much more dramatic for the deflagration speed. In 
particular, for a field strength of i? = 3 T we find the 
deflagration speed Uf = 12.2 m/s for for the magnetic 
field aligned along the easy axis {9 = 0), a much smaller 
speed Uf = 2.6 m/s for 9 = 7r/4 and a negligible value 
Uf = 0.27 m/s for the magnetic field perpendicular to 
the easy axis with 6 = tt/2. Thus we obtain a mag- 
netic deflagration speed almost two orders of magnitude 
smaller for the magnetic field directed along the hard 
axis in comparison with that directed along the easy 
axis. Here we stress that the difference in the deflagra- 
tion speed in our study comes only from modifications 
in the activation energy and Zeeman energy while the 
thermal conduction coefficient remains the same. This 
is different from the experimental studies of Ref. |30 for 
Gd5Ge4 where the deflagration speed changes both be- 
cause of misalignment of the magnetic fleld and thermal 
conduction simultaneously. As a result, the geometry 
suggested here provides better conditions for investigat- 
ing quantum-mechanical properties of the nanomagnets 
(i.e., magnetic anisotropy) and thermal properties of the 
crystals separately. We also stress that the present nu- 
merical results rely on the available models for the nano- 
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FIG. 13. (Color online) Stationary profiles of the scaled tem- 
perature r, fraction of molecules in the metastable state n, 
pressure P (for detonation), and scaled energy release for (a) 
deflagration and (b) detonation for // = 3 T. The charac- 
teristic length scales are Lj = n/Uf — 1.4 /xm for magnetic 
deflagration and Lo = cqtr ~ 0.2 mm for magnetic detona- 
tion. 



magnet Hamiltonian for Mni2-acetateji^ by modifying 
the coefRcients in the Hamiltonian one comes to other nu- 
merical values for the magnetic deflagration speed. The 
present work may also serve for solving the inverse prob- 
lem: by comparing the numerical predictions to future 
refined experiments one may adjust the coefficients in 
the Hamiltonian for nanomagnets. 



B. Magnetic detonation 

The same method may also be used to investigate 
anisotropic properties of magnetic detonation. In con- 
trast to deflagration, magnetic detonation propagates 
due to a leading shock wave preheating the initially 
cold crystal, see Fig. [T3] for typical profiles of tempera- 
ture, pressure and fraction of molecules in the metastable 
state. For comparison, in magnetic deflagration, pre- 
heating happens due to thermal conduction, which is 
negligible for the fast process of magnetic detonation. 
Another important feature of Fig. [T3] (a) is that the 



preheating zone for magnetic deflagration is compara- 
ble by width to the zone of spin switching and energy 
release at H = 3T. This is qualitatively different from 
the analytical Zeldovich-Frank-Kamenetsky deflagration 
model, ^^'■^'^ which assumes a wide preheating region and 
an extremely narrow zone of energy release. We also 
point out that magnetic detonation is noticeably differ- 
ent from the common detonation model (the Zeldovich- 
von Neumann-Doring model) employed in combustion 
science. In particular, the combustion model involves 
a strong delay of the energy release behind the lead- 
ing shockj^ In contrast to that, in magnetic detona- 
tion the spin switching and energy release start directly 
at the leading shock at = 3 T. The most impor- 
tant properties of magnetic detonation propagating along 
the easy axis have been studied in Ref. 22. In par- 
ticular, Ref. [12 has demonstrated that magnetic deto- 
nation is ultimately weak in comparison with common 
combustion detonations21 and, therefore, it propagates 
with a velocity only slightly exceeding the sound speed 
(co ~ 2000 m/s for Mni2-acetate). As a result, the mag- 
netic detonation speed does not depend on the direction 
of the magnetic field. Unlike that, other properties of 
magnetic detonation are quite sensitive to the energy re- 
lease in the spin switching and hence to the magnetic 
field direction. This dependence concerns first of all the 
temperature behind the leading shock (label s), which 
may be calculated as2^ 



T^+' = {a + l)im+ 1) 



3AkBCo 



2TQ 

m + I 



3/2 



(45) 



where F « 2 is the Gruneisen coefficient, and the factor 
m « 4 characterizes the elastic contribution to the pres- 
sure P oc (p/po)™ — 1, where po « 1.38 x lO'^kg/m'^ is 
the initial density of the crystal, see Ref. for details. 
The temperature behind the magnetic detonation front 
(labeled d) depends also on the Zeeman energy release as 



Q 



1 



12co 



2rQ_ 

m + I 



3/2' 



(46) 



The characteristic times of spin switching in magnetic 
detonation at the shock, Tg ~ Ti? exp (Sa/Tg), and at the 
final detonation temperature, ~ Tjiexp{£a/Td), are 
also strongly influenced by the direction of the magnetic 
field. The anisotropic dependence of the temperature on 
the angle between the magnetic fleld and the easy axis 
is presented in Fig. [TH Similarly to deflagration, the 
temperature in magnetic detonation exhibits noticeable, 
though not dramatic, decrease with the angle between 
the magnetic field and the easy axis. For example, for 
H — 5T, the temperature just behind the shock changes 
from Ts = 6.38 K at 6* = to = 4.9 K at 6* = 7r/2; the 
resulting temperature behind the magnetic detonation 
front changes from Ta = 18.7 K at 6^ = to = 15.7 K 
at 9 = 7r/2. However, these moderate modifications of 
temperature, together with respective modifications of 
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the activation energy, lead to dramatic changes in the 
characteristic spin-switching time at the shock, Tg, and 
at the final detonation temperature, r^. For example, for 
the same magnetic field strength as used in the above 
example, H — 5T, we find the reversal time behind the 
leading shock Ts — 2.2 x 10^^ s at 6* = and Tg = 2.3 s 
at 9 = 7r/2; thus we observe variations of the reversal 
time by six orders of magnitude. Such an increase of the 
spin-reversal time makes the magnetic detonation front 
unrealistically wide at large angles so that magnetic det- 
onation becomes impossible for noticeable misalignment 
between the external magnetic field and the easy axis. 
The characteristic spin-switching time at the final tem- 
perature Td and the external field H = 5T changes from 
Td = 2.8 X 10"^ s at 6* = to Trf = 2 X 10^*^ s at 6* = 7r/2. 
Note that in Fig. [T3] we consider larger values of the ex- 
ternal magnetic field than those used in the magnetic 
deflagration experiments. As pointed out in Ref. [22l . 
moderate magnetic fields lead to a quite large thickness 
of the magnetic detonation front, ~ cqT,, much larger 
than the typical sample size unless the magnetic detona- 
tion is formed at a specific resonant field characterizing 
nanomagnets-^" Investigation of spin avalanches at the 
resonant field requires further work beyond the scope of 
the present paper. 



SUMMARY 



FIG. 14. (Color online) Temperature behind the leading 
shock and behind the detonation front (a) versus the external 
magnetic field, and (b) versus the angle between the magnetic 
field and the easy axis 
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FIG. 15. (Color online) The characteristic spin-switching time 
at the shock and at the final detonation temperature versus 
the angle between the magnetic field and the easy axis. 



In this paper, we have investigated anisotropic prop- 
erties of spin avalanches in crystals of nanomagnets 
propagating in the form of pseudo-combustion fronts - 
magnetic deflagration and detonation. In general, the 
anisotropy is expected to be of two types: magnetic 
and thermal. We have focused here on the magnetic 
anisotropy related to the misalignment of the external 
magnetic field and the easy axis of the crystal. The ther- 
mal anisotropy is not considered since at present there is 
no sufficient experimental data for such a study. 

The magnetic anisotropy affects primarily two values 
of the key importance for the magnetic deflagration and 
detonation dynamics - the activation energy and the Zee- 
man energy. Here, we calculated the activation and Zee- 
man energies as a solution to the quantum-mechanical 
problem of a single nanomagnet of Mni2-acetate placed 
in the external magnetic fleld, which is then reversed. 
We demonstrated strong dependence of the activation 
and Zeeman energies on the strength and direction of 
the external magnetic fleld. 

We obtain thated, because of this strong dependence, 
the magnetic deflagration speed is quite sensitive to the 
direction of the magnetic field too. In particular, we 
found that the magnetic defiagration speed may decrease 
by two orders of magnitude for the magnetic field aligned 
along the hard crystal axis instead of the easy one. 

In contrast to magnetic defiagration, the magnetic det- 
onation speed is determined mostly by the sound speed in 
the crystal and, hence, does not depend on the direction 
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of the magnetic field. At the same time, other properties 
of magnetic detonation, such as the temperature behind 
the leading shock and for completed spin reversal, and 
the characteristic time of spin switching, demonstrate a 
strong anisotropy. 
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